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Abstract 

In this article, the nonlinear damping characteristics of magnetorheological damper are expressed with hyperbolic tan¬ 
gent model to simulate its mechanical experimental results. The fitted hyperbolic tangent model can represent hysteretic 
behavior for magnetorheological damper exactly. Based on the hyperbolic tangent model, a quarter-car model with mag¬ 
netorheological damper is established, and a new hybrid fuzzy and fuzzy proportional-integral-derivative (HFFPID) con¬ 
troller integrated with hybrid fuzzy control and fuzzy proportional-integral-derivative control is developed to improve 
the semi-active suspension performance, which can overcome the absence of precise mathematical model. Furthermore, 
numerical simulations for fuzzy proportional-integral-derivative (PID), hybrid fuzzy proportional-integral-derivative 
(HFPID), and HFFPID controllers are investigated to demonstrate the effectiveness of the proposed approaches. The simu¬ 
lation results show that the body acceleration, suspension deflection, and tyre displacement can be reduced more effectively 
using HFFPID controller under sinusoidal road excitation. It can be further concluded that the suspension performance is 
improved more effectively by using HFFPID controller under random road excitation, especially in the peak points. 
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Introduction 

Suspension is one of the most important parts for vehi¬ 
cle, which directly affects the vehicle performance such 
as ride comfort and handling stability. To improve the 
damping effect of vehicle suspension, the traditional 
passive suspension, active suspension, and semi-active 
suspension appeared in succession . 1-5 Traditional pas¬ 
sive suspension is composed of elastic elements and 
damping elements with fixed parameters, which has 
been widely used in vehicle suspension since it requires 
no external energy input and simple structure. 
However, it is only an optimal compromise scheme 
with poor adaptabilities to the varied driving condi¬ 
tions and arbitrary road excitations . 6 Active suspension 
can obtain a high-quality vibration isolation system to 


achieve the ideal control for vehicle through a control¬ 
lable actuator, which has great potential in improving 
the suspension performance. However, its development 
and apphcation are constrained due to its large energy 
consumption, high cost, and complex structure . 7 Semi¬ 
active suspension can accommodate to various roads 
and driving conditions to improve ride comfort and 
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handling stability by changing the damping characteris¬ 
tics of the shock absorber, which is close to active sus¬ 
pension in the vibration attenuation effect. It has 
attracted more attention in numerous scholars and 
manufacturers because of its advantages including high 
energy efficiency and low cost . 8 

Magnetorheological (MR) damper is an intelligent 
device extensively applied in semi-active suspension sys¬ 
tems to damp the undesired vibrations. There is a mag¬ 
netic field generated by applying a direct current to the 
coil in MR damper, which makes MR fluids changing 
from the viscous liquid state to the semi-solid state in 
the resistance gap, thereby a controllable damping force 
can be obtained. MR damper has many distinct advan¬ 
tages such as simple structure, small size, low energy 
consumption, fast response, large dynamic range, and 
continuously controllable resistance . 9 However, the 
damping characteristics of MR damper exhibit strong 
nonlinear hysteretic characteristics. Therefore, it is one 
of the key tasks to establish a more accurate mechanical 
model of MR damper for the design of semi-active sus¬ 
pension system. To describe the dynamic damping 
characteristics accurately, many experts proposed dif¬ 
ferent kinds of mechanical models for MR damper, 
which can be divided into non-parametric and para¬ 
metric models. Among these, non-parametric models 
mainly include polynomial model, neural network 
model, and adaptive neuro fuzzy model . 10-12 These 
models can trace the dynamic behavior for MR damper 
without considering the general physical characteristics 
of implemented parameters. However, parametric mod¬ 
els need a physical modeling process, which is com¬ 
posed of certain mechanical elements including 
dashpot, spring, and friction. Parametric models 
applied in the semi-active control system mainly include 
Bingham model. Phenomenological model, Bouc-Wen 
model, hyperbolic tangent model, and Dahl model . 13-17 
In recent study, the hyperbolic tangent model is selected 
to fit the experimental results of MR damper. This 
selected model can also be effectively used for designing 
inverse model for MR damper. 

Another key task is to design the corresponding con¬ 
trol algorithm for semi-active suspension system. An 
appropriate control algorithm can make semi-active 
suspension to achieve the same vibration attenuation 
effect as active suspension. However, heavy nonlinear 
behaviors and modeling uncertainties complicate the 
development of high-performance controllers . 18,19 To 
improve the dynamics of nonlinear and uncertain sus¬ 
pension systems, many advanced nonlinear control 
strategies have been investigated, such as, Piotr Krauze 
et al. investigated the effectiveness of FxLMS algorithm 
in a quarter-car model including the Bouc-Wen model 
of MR damper. It was stated that the FxLMS-based 
algorithm efficiently mitigated vibrations in the range 
of both the low and high frequencies compared with 


the classical skyhook algorithm . 20 Weichao Sun et al. 
addressed the problem of H 00 control for active vehicle 
suspension systems in finite frequency domain. By the 
Generalized Kalman-Yakubovich-Popov lemma, the 
ride comfort has been improved by minimizing the H x 
norm in specific frequency band . 21 H Du et al . 22 pre¬ 
sented a direct voltage control method based on 
Takagi-Sugeno fuzzy algorithm for MR damper. The 
simulation results showed that a continuous control sig¬ 
nal can be obtained and fast control response is guaran¬ 
teed. Balamurugan et al . 23 proposed a generalized MR 
damper model and verified it in vehicle suspension sys¬ 
tem using skyhook sliding mode control algorithm. 
Recently, Devdutt Singh et al. integrated a hybrid fuzzy 
proportional-integral-derivative (HFPID) controller 
with coupled rules controller in combination with MR 
damper in primary suspension system of traveling semi¬ 
active quarter-car system, and polynomial model was 
selected for simulation purpose of experimental results 
of MR damper. Simulation responses showed that the 
proposed controller provided best performance 24 At 
present, many control strategies have been developed 
and evaluated for their practical applicability in semi¬ 
active suspension system, which includes Linear- 
Quadratic-Gaussian (LQG) control, adaptive neuro 
fuzzy control, skyhook on-off control, sliding mode 
control, adaptive robust control, and integrated adap¬ 
tive robust control . 25-35 

HFPID controller is composed of fuzzy controller 
and traditional proportional-integral-derivative (PID) 
controller , 2 which integrates their advantages. 
However, it is difficult to get the precise parameters of 
traditional PID controller in the HFPID controller. 
Thus, the control accuracy is not ensured. To obtain 
optimal performance with strong robustness, a new 
HFFPID controller integrated with hybrid fuzzy con¬ 
trol and fuzzy-PID control is developed on the basis of 
HFPID controller, which can automatically realize 
the optimal tuning of PID parameters based on fuzzy 
inference. Consequently, its adaptive robustness is 
stronger. 

The main objective of this work is related to vibra¬ 
tion control of MR semi-active suspension system with 
different controllers. First, the mechanical experimental 
results of MR damper are modeled using hyperbolic 
tangent model. Then, a quarter-car model with MR 
damper is established. To further validate the effective¬ 
ness of the proposed HFFPID controller, simulation 
results of semi-active suspension system based on 
hyperbolic tangent model controlled by the HFFPID 
controller, HFPID controller, fuzzy-PID controller, 
and passive suspension under sinusoidal road excita¬ 
tion and random road excitation are compared and 
analyzed considering the body acceleration, suspension 
deflection, and tyre displacement as the evaluation 
indices of the vibration attenuation effect. 
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INSTRON machine 


Power supply 


Figure I. Test system of mechanical performance. 



Figure 2. Prototyping of MR damper. 


Hyperbolic tangent model of MR damper 

The MR damper presented in this article is manufac¬ 
tured by our research group, which belongs to shear- 
valve mode damper with single piston rod. Figure 1 
shows the mechanical performance test system for the 
proposed MR damper, and the prototyping of MR 
damper is shown in Figure 2. The test system mainly 
consists of INSTRON tensile machine, WYK-301 type 
DC regulated power supply, control monitor, and the 
proposed MR damper. The upper end of MR damper 
is fixed on the upper beam of test bench, which remains 
motionless. Meanwhile, the lower end of MR damper is 
excited by simulated vibration. The control monitor is 
used to display and save experimental data of the 
damping force. Here, MR damper works in a 1 Hz 
sinusoid excitation with amplitude of 10 mm, and the 
supplied current to MR damper varied from minimum 
0 A to maximum 1A with the step of 0.25 A. The test 
results of force-displacement and force-velocity of MR 
damper are shown in Figure 3. 

The Bingham model has a clear physical meaning 
with simple expression, and only three parameters need 
to be determined by experimental data, but its precision 
is relatively low. Phenomenological model can reflect 
the nonlinear hysteretic characteristics of MR damper 
perfectly, 14 but more parameters need to be identified, 



Figure 3. Damping performances of proposed MR damper: 
(a) force-displacement and (b) force-velocity. 


which is extremely sophisticated. Hyperbolic tangent 
model presented by Kwok et al. 36 consists of the viscous 
damping (dashpot) placed in parallel with spring stiff¬ 
ness and a hysteretic component. The model is given by 

/ = cx + kx + atanhO&t + Ssgn(x)) +/ 0 (1) 

where x and x are the displacement and velocity of the 
piston of MR damper, respectively; c and k are the vis¬ 
cous and stiffness coefficients, respectively; a is a scale 
factor of the hysteresis; /3 and 8 are the parameters con¬ 
trolling the shape of the hysteresis loops; and / 0 is the 
damper force offset. 

Note that the model contains only a simple hyper¬ 
bolic tangent function, so it can be computed efficiently 
by parameter identifications, as well as included and 
implemented in controller design. Therefore, hyperbolic 
tangent model is used to represent the hysteretic beha¬ 
vior of the MR damper. However, the model contains 
a sign function and six parameters needed to be 
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Table 1. Parameter identification results. 

Current (A) 

Parameter 





a i 

a 2 

k 

03 

fo 

0 

39.0427 

753.3670 

0.3614 

160.721 1 

6.9705 

0.25 

74.2446 

735.0597 

0.3704 

304.521 1 

8.0892 

0.5 

97.8829 

826.2937 

0.3296 

407.0751 

5.8681 

0.75 

123.4396 

823.2025 

0.3307 

513.1546 

6.7852 

1 

161.0215 

744.9815 

0.3658 

661.8812 

8.1226 



Figure 4. Relationship between a„ /= I, 3, and current: 
(a) parameter a t and (b) parameter a 3 . 


identified. Furthermore, it is improved to reduce com¬ 
plexity in the model. The identified parameters are 
reduced to five, and the governing equation for this 
improved model is given as 

/ = «i tanh (a 2 (x + he)) + a 3 (x + kx) +/ 0 (2) 


where a x is a scale factor of the hysteresis, which is 
related to the control current; a 2 and a 3 are the para¬ 
meters related to the viscous damping coefficient in the 
pre-yield region to post-yield region, respectively; and 
k = V 0 /X 0 (V 0 and X 0 are defined as the absolute value 
of critical piston velocity and critical piston displace¬ 
ment when the damping force / is zero, respectively), 
which is a scale factor of hysteresis loop width. 

It appears that the maximum damping force of MR 
damper increases gradually with the increasing control 
current as shown in Figure 3. The coefficients in equa¬ 
tion (2) under the different current are obtained 
through parameter identification solved by least-square 
method based on the above experimental results of MR 
damper. The principle of identification is to acquire the 
unknown five coefficients to minimize equation (3), 
that is, to minimize the error sum of squares between 
the predicted force and the measured force 

e = El = 1 I fexperiment ~ fpredict] (3) 

Parameter identification results are given in Table 1. 
Observing Table 1, the coefficients associated with cur¬ 
rent in equation (2) can be distinguished by studying 
the relationship between the same coefficients and cur¬ 
rent under the different currents. 

Table 1 shows that the parameters a 2 , k, and / 0 are 
not changed significantly under the different current, 
which can be considered as a constant value, and their 
mean value a 2 = 776.5809, k = 0.3506, and/ 0 = 7.1671 
are selected. Parameters a\ and a 3 increase along with 
the increasing current, which have an evident linear 
relationship. The relationship between these two vari¬ 
ables and current is shown in Figure 4. 

The relationship between parameters and a 3 , and 
current are obtained by linear regression approach, and 
it is expressed as 

( a x = b x I + ci = 117.2610/ + 40.4957 . , 

\ a 3 = b 2 I + c 2 = 484.3815/ + 167.2799 * ’ 

Therefore, hyperbolic tangent model is given by 
/ = (b x I + ci) tanh (a 2 (x + kx)) + (b 2 I + c 2 )(x + he) + f 0 

( 5 ) 
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Figure 5. Comparison of predicted and experimental results: 
(a) force-displacement and (b) force-velocity. 

A comparison between the predicted responses for 
hyperbolic tangent model and experimental results 
under different currents are provided in Figure 5. It 
can be concluded that the improved hyperbolic tangent 
model matches very closely with the experimental 
results, which can be effectively applied to semi-active 
control. Moreover, the desired damping force supplied 
by controller as well as the piston displacement and 
velocity can be known easily. 



Table 2. Parameter values of the quarter-car suspension 
model. 


Parameter 

m b (kg) 

(kg) 

k s (N/m) 

k, (N/m) 

Value 

219 

53 

12,480 

240,000 


The governing equations of motion for the sprung mass 
and unsprung mass are 


/ rribz b + k s (z b - z w ) -/ = 0 
\ m w z w + k s (z w - z b ) + k t (z w - z r ) +f = 0 


( 6 ) 


where ny, is the sprung mass, which represents the car 
body; m w is the unsprung mass, which represents the 
wheel assembly; z b and z w are the displacements of the 
sprung mass and unsprung mass, respectively; z r is the 
road displacement input; k s and k, are the suspension 
stiffness and tyre stiffness, respectively; and / represents 
the damping control force generated by MR damper. 
The parameter values of the quarter-car suspension 
model are listed in Table 2. 

Taking dynamic relationship into account, the road 
displacement z r and damping control force/ are selected 
as the input vectors, and the state variables can be 
defined as 


Quarter-car suspension model 

In this section, a quarter-car suspension model with 
MR damper is taken for comparative study, as shown 
in Figure 6. 

This model is selected to investigate vertical 
dynamics of the complete structure with neglecting 
other motions such as vehicle roll and pitch dynamics. 


^ ^ (?) 

Accordingly, a state-space equation for equation (6) 
is derived as 

z = Az + Bu (8) 

where 
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Table 3. Rule base for fuzzy controller. 

dele 

NB 

NM 

NS 

ZE 

PS 

PM 

PB 

NB 

NM 

NM 

NM 

NM 

NM 

NS 

ZE 

NM 

NM 

NM 

NM 

NM 

NS 

ZE 

PS 

NS 

NM 

NM 

NM 

NS 

ZE 

PS 

PM 

ZE 

NM 

NM 

NS 

ZE 

PS 

PM 

PM 

PS 

NM 

NS 

ZE 

PS 

PM 

PM 

PM 

PM 

NS 

ZE 

PS 

PM 

PM 

PM 

PM 

PB 

ZE 

PS 

PM 

PB 

PM 

PM 

PM 


NB: negative big; NM: negative medium; NS: negative small, ZE: zero; PS: positive small; PM: positive medium; PB: positive big. 



Semi-active control algorithm 

Semi-active controller design 

First, the fuzzy controller does not require the accurate 
mathematical model for the controlled object, and it 
can achieve the effective nonlinear control effect 
according to the implemented operator’s experience 
and knowledge, which has been frequently tested in 
various active and semi-active suspensions for vibration 
control. 37,38 Second, the fuzzy-PID controller was often 
used in vehicle suspension vibration control due to its 
strong adaptability and robust performance. 39 The 
F1FFPID controller combines the advantages of both 
fuzzy controller and fuzzy-PID controller in a single 
hybrid-type controller. 

Fuzzy controller. Fuzzy controller has two input and one 
output variables. One of the input variables chosen for 
the fuzzy controller is the error (e) between the car 
body velocity and its reference value, and the change 
rate of the error ( de ) denoted by the error between the 
car body acceleration and its reference value is selected 
as another input variable. The desired damping force 
(w) is selected as the output variable. The seven linguis¬ 
tic variables that are selected to describe the input and 
output variables are as follows: NB (negative big), NM 
(negative medium), NS (negative small), ZE (zero), PS 
(positive small), PM (positive medium), and PB 



Figure 8. Fuzzy controller scheme. 


(positive big). In this article, triangular membership 
functions (MFs) are used for linguistic variables. Thus, 
the concerned MFs for the input and output variables 
are shown in Figure 7. Variables are all defined in the 
interval [—1, 1]. Actual intervals of variables are 
decided by the following scaling factors: k v , k a , and k u , 
respectively. Fuzzy controller scheme is shown in 
Figure 8. 

According to the fuzzy rules of the rule base look-up 
table: if e is E t and de is Ej, then u is Vyj), the corre¬ 
sponding fuzzy control rule base is presented in 
Table 3. For instance, one of the possible rules is that if 
e is NS and de is PB, then u is PM. This rule can be 








































Hu et al. 



Figure 9. Fuzzy-PID controller scheme. 



Figure 10. MFs of the output variables. 


explained as follows: if the error is small, the car body 
velocity is around the reference value. Meanwhile, a 
positive big value of change rate of the error shows that 
the car body velocity rapidly approaches to the refer¬ 
ence value, but the car body acceleration is partial 
large. Consequently, the controller output should be a 
positive medium value to prevent overshoot and make 
a brake effect. Moreover, Mamdani method is selected 
as the fuzzy inference method, and the centroid method 
is used for defuzzification, which converts the linguistic 
variables to numerical data. 

Fuzzy-PID controller. In the conventional PID controller, 
the ideal characteristics of the proportion, integration, 
and differentiation can be expressed as 

u(t) = K p e(t) + K t I e(t)dt + K d d -§ (9) 

o 

where K p is the proportional gain; K t is the integral 
gain; K d is the differential gain; e(t) represents the con¬ 
troller input, which is the error between the car body 
velocity and its reference value; and u(t) represents the 
controller output, which is the desired damping force. 

Fuzzy-PID controller is a composite controller 
which connects with fuzzy control and conventional 
PID control, as shown in Figure 9. It can automatically 
modify the three parameters online for PID controller 
according to the fuzzy control theory to meet the differ¬ 
ent requirements, which merges the advantages of fuzzy 
controller and PID controller including flexibility, 
robustness, simplicity, and reliability. So, this strategy 
improves the dynamic response of the system with high 
dynamic and static performance. 40 Observing Figure 9, 
the specific working process of the fuzzy-PID controller 
can he explained as follows: first, the output variables 
A K p , A K h and A K d are obtained by fuzzy controller 
with e and de as the input variables through fuzzy infer¬ 
ence. Subsequently, they are all output to the PID con¬ 
troller so that the PID parameters are tuned 
dynamically in real time. Finally, the PID control para¬ 
meters K p , K h and K d are received. It can be repre¬ 
sented as 


K p = Kpo + A K p Xq p 
Ki = K i0 + AKi X q t (10) 

K d = K m + A K d X q d 

where K p , K b and K d are the final parameters of fuzzy- 
PID controller; K p0 , K i0 , and K d0 are the initial para¬ 
meters of PID controller; A K p , A K b and A K d are the 
output variables of fuzzy controller; and q p , </,, and q d 
are the corresponding correction coefficient. 

The above seven linguistic variables are used to 
describe the input variables in this fuzzy controller, and 
the MFs are also shown in Figure 7. The output vari¬ 
ables are all portrayed by four linguistic variables: ZE 
(zero), S (small), M (medium), and B (big), as shown in 
Figure 10. The output variables are all defined in inter¬ 
val [0, 1]. According to the rule base look-up table, if e 
is E{ and de is Ej, then u is Uqj), the corresponding 
fuzzy control rule bases are presented in Tables 4-6, 
separately. Equally, Mamdani method is selected as the 
fuzzy inference method, and the centroid method is 
used for defuzzification in this fuzzy controller. 

HFFPID controller. An intelligent switching mechanism is 
added for HFFPID controller. Conventional PID con¬ 
troller in HFPID controller is replaced by fuzzy-PID 
controller with stronger adaptive robustness in 
HFFPID controller. It can be automatically switched 
between fuzzy controller and fuzzy-PID controller 
according to whether the error is larger than the thresh¬ 
old. Its detailed principle can be described as follows 

\,if\e\>w Fuzzy controller 
0, if\e\ <w Fuzzy-PID controller ^ ^ 

where w is the switching threshold value. 

The fuzzy controller is chosen as the output of the 
control system when the car body velocity is larger than 
the switching threshold value. In this case, fuzzy con¬ 
troller has better robustness. Otherwise, if the car body 
velocity is smaller than the switching threshold value, 
the control system transfers immediately the controller 
from fuzzy controller to fuzzy-PID controller. In 
the car body velocity near the reference value, the 
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Table 4. 

Rule base for variable AKp. 







dele 

NB 

NM 

NS 

ZE 

PS 

PM 

PB 

NB 

B 

B 

M 

M 

S 

S 

ZE 

NM 

B 

M 

M 

S 

S 

M 

ZE 

NS 

B 

M 

S 

ZE 

S 

S 

B 

ZE 

B 

S 

ZE 

ZE 

ZE 

s 

B 

PS 

B 

S 

S 

ZE 

S 

M 

B 

PM 

ZE 

ZE 

S 

S 

M 

M 

B 

PB 

ZE 

ZE 

S 

M 

M 

B 

B 

NB: negative big; NM: negative medium; NS: negative small, ZE: 

Table 5. Rule base for variable AK,. 


; PS: positive small; PM: positive medium; 

PB: positive big. 


dele 

NB 

NM 

NS 

ZE 

PS 

PM 

PB 

NB 

ZE 

S 

M 

B 

M 

S 

ZE 

NM 

ZE 

S 

B 

B 

B 

S 

ZE 

NS 

ZE 

ZE 

M 

B 

M 

ZE 

ZE 

ZE 

ZE 

ZE 

B 

B 

B 

ZE 

ZE 

PS 

ZE 

ZE 

B 

B 

B 

ZE 

ZE 

PM 

ZE 

S 

M 

B 

M 

S 

ZE 

PB 

ZE 

S 

M 

B 

M 

S 

ZE 

NB: negative big; NM: negative medium; NS: negative small, ZE: 

Table 6. Rule base for variable A K d . 


; PS: positive small; PM: positive medium; 

PB: positive big. 


dele 

NB 

NM 

NS 

ZE 

PS 

PM 

PB 

NB 

S 

M 

ZE 

ZE 

ZE 

M 

S 

NM 

S 

M 

S 

ZE 

S 

M 

S 

NS 

S 

B 

S 

S 

S 

B 

S 

ZE 

S 

B 

S 

ZE 

S 

B 

S 

PS 

S 

B 

S 

S 

S 

B 

S 

PM 

S 

M 

S 

ZE 

S 

M 

S 

PB 

S 

M 

ZE 

ZE 

ZE 

M 

S 


NB: negative big; NM: negative medium; NS: negative small, ZE: zero; PS: positive small; PM: positive medium; PB: positive big. 


fuzzy-PID controller has better accuracy and effective¬ 
ness. Thus, the designed HFFPID controller takes full 
advantage of both the fuzzy controller and fuzzy-PID 
controller. 

Semi-active control constraint 

The damping force generated by MR damper under 
the same piston velocity and displacement is a bounded 
value because the current inputs into the MR damper 
coil is between the minimum current and the maximum 
current. Accordingly, it is necessary to saturate the 
desired damping control force obtained by the semi¬ 
active controller. The constrained damping control 
force is expressed as 

|/ m i„W|<UimW|<|/maxW| (12) 


where f m i n (t) is the minimum damping force at a certain 
piston velocity; f max (t) is the maximum damping force 
at a certain piston velocity; and /ii m (t) is the constrained 
control damping force. The maximum input current of 
MR damper is 1 A, and the minimum damping force 
and the maximum damping force can be calculated at a 
certain piston velocity and displacement according to 
formula (5), respectively. 

It can be seen from Figure 5(b) that the damping 
force provided by MR damper is basically located in 
the first and third quadrants in the force-velocity curve. 
Hence, MR damper can output the damping force to 
track the desired control force when the symbol of the 
desired control force is the same as that of the piston 
velocity. However, MR damper can hardly track the 
desired control force when the desired control force 
and piston velocity are opposite sign, and the input 
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Figure I I. Control schematic for semi-active suspension 
system. 



Figure 12. Body acceleration under sinusoidal road excitation. 


current is set to 0 A under this circumstances. The final 
traceable damping control force can be obtained as 


ftrait) = { 


/limW 

/min (0 


fe{t)-v{t)< 0 


(13) 


where f tr Jf) represents the traceable damping control 
force; f e {t) represents the desired control force; and v(t) 
is the piston velocity. 

According to formulas (5), (12), and (13), the corre¬ 
sponding control current in the hyperbolic tangent 
inverse model is 


= ftrait) - c\m - c 2 n -f 0 

b\m + & 2 « 


(14) 


where m = tanh(a 2 (x(t) + kx{t))), n = x(t) + kx{t). 



Numerical simulation analysis 

To validate the effectiveness of the designed HFFPID 
controller and evaluate various suspension perfor¬ 
mance, two types of road excitations are considered in 
the simulation. Considering the car body acceleration, 
suspension deflection, and tyre displacement as the eva¬ 
luation indices of the vibration attenuation effect, simu¬ 
lation of semi-active suspension is conducted by the 
hyperbolic tangent model. Comparisons among differ¬ 
ent controllers, containing the HFFPID controller, 
HFPID controller, fuzzy-PID controller, and passive 
suspension (i.e. input control current is OA) are ana¬ 
lyzed. The parameter values of the quarter-car suspen¬ 
sion model are listed in Table 2. Control schematic for 
semi-active suspension system is shown in Figure 11. 

A sinusoidal road input is adopted as the first road 
excitation with amplitude of 0.07 m and a frequency of 
1 Hz. The body acceleration A b , suspension deflection 
D s , and tyre displacement D, under sinusoidal road 
excitation are shown in Figures 12-14, respectively. 


Figure 13. Suspension deflection under sinusoidal road 
excitation. 


V 


V 


Figure 14. Tyre displacement under sinusoidal road excitation. 
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Figure 15. Comparison of the actual damping force and the 
desired control force under sinusoidal road excitation: (a) fuzzy- 
PID, (b) HFPID, and (c) HFFPID. 
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Figure 16. Input control current under sinusoidal road 
excitation: (a) fuzzy-PID, (b) HFPID, and (c) HFFPID. 


It can be seen that the better vibration reduction per¬ 
formance are obtained for the controlled semi-active 
suspensions compared to the uncontrolled passive sus¬ 
pension, and HFFPID controller has the best vibration 
attenuation performance compared with other control¬ 
lers in semi-active suspensions. In semi-active suspen¬ 
sion control, the comparison between the actual 
damping control force generated by MR damper and 


the desired control force generated by controllers is 
provided in Figure 15, and the input control current 
for MR damper under different controllers is provided 
in Figure 16. It is also observed that the desired control 
force is relatively smooth, whereas the mutation phe¬ 
nomenon appears in the actual damping control force 
and input control current. The reason for this is that 
the input control current is suddenly loaded when the 
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Figure 17. Reduction percentage of RMS under sinusoidal 
road excitation. 


desired control force is satisfied from the non- 
constrained state to constrained state. But it can also 
be concluded from Figure 15 that both of them show 
the extremely similar variation trend, which means that 
the hyperbolic tangent model of MR damper can pre¬ 
cisely track the desired control force so that the semi¬ 
active suspension performance is almost equal to the 
active one. According to equation (14), the input con¬ 
trol current controlled by fuzzy-PID controller is rela¬ 
tively small because its suspension deflection and 
relative velocity are relatively large. Conversely, the 
input control current controlled by HFFPID controller 
is relatively large because its suspension deflection and 
relative velocity are smaller than other controllers. 

Compared with passive suspension, the reduction 
percentage of root mean square (RMS) value of the 
body acceleration, suspension deflection, and tyre 
displacement under sinusoidal road excitation for 
fuzzy-PID controller, HFPID controller, and HFFPID 
controller is shown in Figure 17. The body acceleration 
controlled by HFFPID controller, HFPID controller, 
and fuzzy-PID controller is reduced by 58.9%, 42.1%, 
and 37.3% compared with the passive suspension, 
respectively, and their suspension deflection are reduced 
by 58.7%, 40.2%, and 42.2%, respectively. Their tyre 
displacements are reduced by 22.7%, 17.8%, and 7.8%, 
respectively. It can be concluded that the body accelera¬ 
tion, suspension deflection, and tyre displacement can 
be reduced more effectively by HFFPID controller 
compared with passive suspension and other control¬ 
lers, which results in a better vibration attenuation 
effect. In addition, the body acceleration and tyre dis¬ 
placement controlled by HFPID controller are smaller 
than fuzzy-PID controller, while the suspension deflec¬ 
tion performance controlled by fuzzy-PID controller is 
better than HFPID controller. 



0 1 2 3 4 5 

Time (s) 


Figure 18 . Body acceleration under random road excitation. 


To further validate the effectiveness of the designed 
HFFPID controller, a second type of road excitation, 
that is, a random road excitation, is applied. 41 The road 
displacement of the random road excitation is gener¬ 
ated by 


z r {t) = - a ■ v • z r (t) + 8 ■ ■ wo(t) (15) 

where z r (t ) is the road displacement; a and 8 are the 
undetermined coefficients; G g (n 0 ) is the road irregular¬ 
ity coefficient; w o (0 is the white noise with covariance 
of 1; and v is the vehicle speed. 

In the simulation, a = 0.111m -1 , 8 = 4.44 m, 
G q (n 0 ) = 1024 X 10 _6 m 3 in D grade road, and 
v = 20 m/s. The body acceleration and suspension 
deflection under random road excitation are presented 
in Figures 18 and 19, respectively. Further analysis 
reveals that the proposed HFFPID controller is effec¬ 
tive to improve suspension performance. Meanwhile, 
the comparison between the actual damping control 
force and desired control force is provided in 
Figure 20, and the input control current for MR dam¬ 
per is provided in Figure 21. The simulation results 
for this case prove the validity of the hyperbolic tan¬ 
gent inverse model again. 

Under random road excitation, the reduction per¬ 
centage of the RMS value of the body acceleration and 
suspension deflection compared with passive suspen¬ 
sion is presented in Figure 22, respectively, for fuzzy- 
PID controller, HFPID controller, and HFFPID 
controller. Figure 22 shows that the body acceleration 
controlled by HFFPID controller, HFPID controller, 
and fuzzy-PID controller are decreased, respectively, 
by 55.9%, 42.6%, and 37.9% compared with the pas¬ 
sive suspension, and their suspension deflection is 
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excitation. 


decreased, respectively, by 49.8%, 22.5%, and 17.5%. 
It can be inferred that the body acceleration and sus¬ 
pension deflection controlled by HFFPID controller 
are less than other controllers, especially in the peak 
points. Moreover, the body acceleration and suspen¬ 
sion deflection controlled by HFFPID controller still 
maintain a similar reduction percentage compared with 
that in sinusoidal road excitation, whereas the reduc¬ 
tion percentage of the suspension deflection controlled 
by fuzzy-PID controller and HFFPID controller are 
decreased remarkably compared with that in sinusoidal 
road excitation. The results show that HFFPID con¬ 
troller has better robustness, which can provide the 
advanced suspension performance and improve ride 
comfort in general. 

Conclusion 

In the semi-active suspension control system using MR 
damper, the vibration attenuation effect depends on 
not only the control algorithm of the system but also 
the accurate mechanical model of MR damper. In this 
article, the mechanical experimental results of MR 
damper are fitted using hyperbolic tangent model, and 
a quarter-car model with MR damper is also estab¬ 
lished. At the same time, a novel HFFPID controller 
integrated with fuzzy-PID control and fuzzy control is 
developed. Furthermore, considering the body accel¬ 
eration, suspension deflection, and tyre displacement 
as the evaluation indices of the vibration attenuation 
effect, the control performances using the proposed 
hybrid controller are simulated compared with 
HFPID controller, fuzzy-PID controller, and passive 
suspension. 




(b) 

1 


— “1 | 



Figure 20. Comparison of the actual damping force and the 
desired control force under random road excitation: (a) fuzzy- 
PID, (b) HFPID, and (c) HFFPID. 

Under sinusoidal road excitation, the results show 
that the body acceleration of semi-active suspension 
controlled by HFFPID, HFPID, and fuzzy-PID con¬ 
troller is reduced by 58.9%, 42.1%, and 37.3%, respec¬ 
tively, compared with the passive suspension, and their 
suspension deflection are reduced by 58.7%, 40.2%, 
and 42.2%, respectively, and their tyre displacements 
are reduced by 22.7%, 17.8%, and 7.8%, respectively. 
The tyre displacement reduction is not obvious due to 
the contradiction between the ride comfort and 
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Figure 21. Input control current under random road 
excitation: (a) fuzzy-PID, (b) HFPID, and (c) HFFPID. 
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Figure 22. Reduction percentage of RMS under random road 
excitation. 


HFPID, and fuzzy-PID controller are decreased by 
55.9%, 42.6%, and 37.9%, respectively, compared with 
the passive suspension, and their suspension deflection 
are decreased by 49.8%, 22.5%, and 17.5%, respec¬ 
tively. It can be further seen that the HFFPID control¬ 
ler improves more effectively the vibration reduction 
performance of suspension than other controllers, espe¬ 
cially at peak points, which gives better ride comfort. 
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handling stability control. Simultaneously, the actual 
output damping force has a similar varied trend with 
the expected control force in the semi-active suspension 
control. Thus, the hyperbolic tangent model of MR 
damper can track the expected control force. 

Under random road excitation, the body accelera¬ 
tion of semi-active suspension controlled by HFFPID, 
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